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One of the long-standing problems in antenna engineering is the realization of highly-directive 
beams using low-profile devices. In this paper we provide a solution to this problem by means of 
Huygens’ metasurfaces (HMSs), based on the equivalence principle. This principle states that a 
given excitation can be transformed to a desirable aperture field by inducing suitable electric and 
magnetic surface currents. Building on this concept, we propose and demonstrate cavity-excited 
HMS antennas, where the single-source cavity excitation is designed to optimize aperture illumina¬ 
tion, while the HMS facilitates the current distribution that ensnres phase parity of aperture fields. 

The HMS breaks the coupling between the excitation and radiation spectrum typical to standard 
partially-reflecting surfaces, allowing tailoring of the aperture properties to produce a desirable ra¬ 
diation pattern. As shown, a single semianalytical formalism can be followed to achieve control of 
a variety of radiation features, such as the direction of the main beam or the side lobe level, by 
proper modification of the HMS and the source position. Relying on a cavity excitation, this can 
be achieved without incurring edge-taper losses and without any degradation of the aperture illu¬ 
mination for arbitrarily-large apertures. With the recent demonstrations of Huygens’ metasurfaces 
at microwave, terahertz, and optical frequencies, the proposed low-profile design may find its use 
in a myriad of applications across the electromagnetic spectrum, from highly-directive antennas to 
highly-efficient quantum-dot emitters, reaching near-unity aperture efficiencies. 

PACS numbers: 41.20.Jb, 78.67.Pt, 84.40.Ba 


I. INTRODUCTION 

Achieving high directivity with compact radiators has 
been a major concern of the antenna community since 
its early days m- Still today, many modern appli¬ 
cations, such as automotive radars, satellite communica¬ 
tion, millimetre-wave point-to-point communication, and 
microwave imaging, strive for simple and efficient low- 
profile antennas producing the narrowest possible beams 
HHi]. Extending the size of the radiating aperture leads 
to an enhanced directivity, but only if the aperture is ef¬ 
ficiently excited. To date, uniform illumination of large 
apertures is achievable with reflectors and lenses; al¬ 
though these can be made compact using concurrent 
metamaterial concepts, they still require substantial sep¬ 
aration between the source and the aperture, resulting 
in a large overall antenna size In addition, feed 

blockage and spillover effects must also be considered, 
usually complicating the design and reducing the device 
efficiency. High aperture efficiencies can also be achieved 
using antenna arrays |10j : nevertheless, the requirement 
for elaborated feed network significantly increases the 
complexity of this solution and limits its compactness 
[ID, and may also introduce considerable feed-network 
losses. 

Leaky-wave antennas (LWAs), on the other hand, can 
produce directive beams using a low-profile structure fed 
by a simple single source |12| . Their typical configu¬ 
ration consists of a guiding structure with a small per¬ 
turbation, facilitating coupling of guided modes to free- 
space radiation. In the much-discussed Fabry-Perot (FP) 


LWAs, a localized source is sandwiched between a per¬ 
fect electric conductor (PEC) and a partially-reflecting 
surface (PRS), forming a longitudinal FP cavity [51 ITl^. 
By tuning the cavity height at the design frequency w, 
favourable coupling of the source to a single parallel-plate 
waveguide mode is achieved, forming a dominant leaky 
wave emanating from the source; the typical device thick¬ 
ness lies around half of a wavelength. The leaky mode 
is characterized by a transverse wavenumber whose real 
part kt corresponds to the waveguide dispersion, and is 
accompanied by a small imaginary part a determined 
by the PRS. This leads to an azimuthally-symmetric di¬ 
rective radiation through the PRS towards the direction 
defined by sindout ~ kt/k, where k = ujy/Jie is the free- 
space wavenumber, with a beamwidth proportional to a. 
Broadside radiation is achieved when 0out is small enough 
such that the splitting condition fct < a is satisfied, and 
the peaks of the conical beam merge m. 

Another class of LWAs which has received significant 
attention lately is based on modulated impedance meta¬ 
surfaces (MoMet As) [TSHTB] . These so-called holographic 
antennas use a point source to excite surface waves on 
a thin dielectric sheet covered with metallic patches and 
backed by a PEC ground plane, establishing effective sur¬ 
face impedance boundary conditions m-, guiding surface 
waves, these structures can be very thin, below fifth of a 
wavelength. Similar to FP-LWAs, the guiding structure 
is designed such that only a single surface mode kf > k 
is allowed to propagate; small modulation of the surface 
impedance, implemented by variation of patch sizes or di¬ 
electric thickness, couples the bound modes to radiative 
modes. To facilitate such coupling in the case of sur- 
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face waves, whose transverse momentum is greater than 
that of free-space, the impedance modulation should have 
a periodicity a = 2 t:/K comparable to the wavelength. 
The interaction between the surface wave and the per¬ 
turbation results in an infinite number of Floquet-Bloch 
(FB) harmonics kt,n = kt+nK] the periodicity should be 
designed such that one of them radiates to the desirable 
direction, while the others become evanescent, ensuring 
good directivity. The leakage rate a, and correspond¬ 
ingly the beamwidth, are determined by the depth of the 
modulation [TS] . 

Both FP-LWAs and MoMetAs have an appealing com¬ 
pact configuration and their radiation characteristics can 
be rather simply controlled by tuning the properties of 
the guiding structure and the perturbation. Nonethe¬ 
less, due to their leaky-wave nature, they suffer from a 
fundamental efficiency limitation when considering prac¬ 
tical finite apertures: designing a moderate leakage rate 
a with respect to the aperture length L yields uniform il¬ 
lumination of the aperture (high aperture efficiency) but 
results in considerable losses from the edges (low radia¬ 
tion efficiency); on the other hand, large values of a lead 
to high radiation efficiencies but in this case only a por¬ 
tion of the aperture is used for radiation, leading to a 
wider beam [2(11 - 122] . 

To mitigate edge-taper losses, shielded FP-LWA struc¬ 
tures have been recently proposed, using PEC side walls 
which form a lateral cavity [2111429] . Nevertheless, the 
tight coupling between the propagation of the leaky mode 
inside the FP cavity and the angular distribution of the 
radiated power manifested by sindout ^ kt/k poses seri¬ 
ous limitations on the achievable aperture efficiency. This 
is most prominent for antennas radiating at broadside, in 
which only low-order lateral modes, carrying transverse 
wavenumbers which are small enough to satisfy the split¬ 
ting condition, can be used. Consequently, such antennas 
are designed to excite exclusively the TEio lateral mode, 
which inherently limits the aperture efficiency, defined as 
the relative directivity with respect to the case of uniform 
illumination, to about 80% [30]. Although this problem 
can be partially solved by the use of artificial magnetic 
conductor (AMC) side walls instead of PECs [3TJ 152] . 
one of the most limiting constraints on the antenna de¬ 
sign which still remains is the requirement for mode pu¬ 
rity. As the dominant spectral components of the field 
in the cavity directly translate to prominent lobes in the 
radiation pattern, only a single mode should be allowed 
to propagate in the cavity to guarantee high directivity. 
However, as demonstrated by [53], suppression of para¬ 
sitic modes in a cavity is usually a very difficult problem. 
In particular, single mode excitation becomes increas¬ 
ingly challenging as the desirable aperture size increases, 
due to the small differences between the wavenumbers 
associated with low-order modes; thus, in practice, this 
solution cannot be used for arbitrarily-large apertures. 

From the discussion so far it follows that it would be 
very beneficial if we could optimize separately the fields 
inside the cavity and the fields formed on the aperture. 


This would allow us to achieve good illumination of the 
aperture without the necessity to meet restricting condi¬ 
tions (e.g., the splitting condition, or single mode exci¬ 
tation), stemming from the coupling between excitation 
and radiation fields. But how to achieve such a sep¬ 
aration? The equivalence principle suggests that for a 
given field exciting a surface, desirable (arbitrary) aper¬ 
ture fields can be formed by inducing suitable electric 
and magnetic surface currents, supporting the required 
field discontinuities [30]. Based on this idea, the concept 
of Huygens’ metasurfaces (HMSs) has been recently pro¬ 
posed as a means for versatile wavefront manipulation 

[5H12|. 

Huygens’ metasurfaces are planar structures composed 
of subwavelength elements (meta-atoms), engineered to 
generate the surface currents required by the equivalence 
principle to achieve a prescribed functionality. In general, 
for a given incident field and desirable transmitted field, 
elements exhibiting effective loss and gain are required 
for the implementation |43L 144] . However, for certain 
applications, the fields can be judiciously stipulated such 
that the metasurface can be constructed from passive and 
lossless elements, i.e. electric and magnetic polarizable 
particles [511 [55]. In fact, we have recently shown that 
if the reflected and transmitted fields are set such that 
the wave impedance and the real power are continuous 
across the two facets of the metasurface, the aperture 
phase can be tailored by a passive and lossless HMS to 
produce directive radiation towards a prescribed angle 
^out j for any given excitation source; the design procedure 
is straightforward once the source plane-wave spectrum 
is assessed [55] . 

Indeed, in this paper we propose to harness the equiv¬ 
alence principle to efficiently convert fields excited in a 
cavity by a localized source to highly-directive radiation 
using a Huygens’ metasurface: cavity-excited HMS an¬ 
tenna. The device structure resembles a typical shielded 
FP-LWA configuration, with an electric line source sur¬ 
rounded by three PEC walls and a Huygens’ metasurface 
replacing the standard PRS (Fig. [^. For a given aper¬ 
ture length L and a desirable transmission angle 0outj we 
optimize the FP cavity thickness and source position to 
predominantly excite the highest-order mode of the lat¬ 
eral cavity, with the HMS reflection coefficient ensuring 
the wave impedance is equalized along the metasurface; 
this guarantees the aperture is well illuminated. Once 
the source configuration is established, we stipulate the 
aperture fields to follow the power profile of the cavity 
mode, ensuring the real power is conserved at each point, 
and impose the suitable linear phase to promote radiation 
towards dout- With the cavity fields and aperture fields 
in hand, we invoke the equivalence principle and evalu¬ 
ate the electric surface impedance and magnetic surface 
admittance required to support the resultant field dis¬ 
continuity [341 [351 113 EZ] • Our previous work [i^ guar¬ 
antees that these would be purely reactive, hence could 
be implemented using passive and lossless meta-atoms. 

Utilizing the equivalence principle as described results 
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v = I/2 


FIG. 1. Physical configuration of a cavity-excited Huygens’ 
metasurface antenna. 


in formation of aperture fields, the magnitude of which 
follows the power distribution inside the cavity, whereas 
their phase is independently determined to vary in a 
plane-wave-like fashion. This has two important implica¬ 
tions. First, as the power profile of the highest-order lat¬ 
eral mode creates hot spots of radiating surface currents 
approximately half a wavelength apart, a uniform virtual 
phased-array is formed on the HMS aperture; based on 
array theory m, such excitation profile is expected to 
yield very high directivity with no grating lobes regard¬ 
less of the scan angle 6>out- Second, in contrast to LWAs 
of any type, the antenna directivity does not deteriorate 
significantly even if other modes are partially excited, as 
these would merely vary the amplitude of the virtual ar¬ 
ray elements, without affecting the phase purity. 

This semianalytical design procedure can be applied to 
arbitrarily-large apertures, yielding near-unity aperture 
efficiencies, in agreement with full-wave simulations; due 
to the PEC side walls, no power is lost via the edges. This 
offers an effective way to overcome the efficiency tradeoff 
inherent to FP-LWAs and MoMetAs, while preserving 
the advantages of a single-feed low-profile antenna. 


II. THEORY 

A. Formulation 


To design the HMS-based antenna, we simply apply 
the general methodology developed in to the source 
configuration of Fig. for completeness, we recall briefly 
its main steps. We consider a 2D scenario {d/dx = 0) 
with the HMS at z = 0 and a given excitation geome¬ 
try at z < z' < 0 emb edded in a homogeneous medium 
(k = ijjy/efl, T] = s/JI/e). Under these circumstances, the 
incident, reflected and transmitted fields in the vicinity of 
the HMS can be expressed via their plane-wave spectrum 

m 

Ar (y, z) = krjIoT-^ {ipfikt) 

< Ef {y, z) = -kyloE-^ {(fc,) / (fc,) (1) 

(y, z) = kyloT-^ {ipT ih) , 


where F ^{g{kt,z)} = ^ dktg {kt] z) is the 
inverse spatial Fourier transform of g(kt\z), / {kt) is the 
source spectrum, F {kt) is the HMS reflection coefficient, 
and T {kt) = T {kt) [1 -I- F {kt)] is the transmission spec¬ 
trum. As before, kt denotes the transverse wavenumber, 
and the longitudinal wavenumber is /3 = yjk'^ — k^. For 
simplicity, we only consider here transverse electric (TE) 
fields {Ez = Ey — = 0); the nonvanishing magnetic 

field components Hy,Hz can be calculated from E^ via 
Maxwell’s equations. 

For a given source spectrum, it is required to determine 
the reflected and transmitted fields, through the respec¬ 
tive degrees of freedom F (fct) and T {kt), that would im¬ 
plement the desirable functionality. Once the tangential 
fields on the two facets of the HMS are set, the equiva¬ 
lence principle is invoked to evaluate the required electric 
and magnetic surface currents to induce them |30j . The 
polarizable particles comprising the HMS are then de¬ 
signed such that the average fields acting on them induce 
these surface currents [461 SI!- Analogously, the HMS 
can be characterized by its electric surface impedance 
Zse {y) and magnetic surface admittance Ysm {y), relat¬ 
ing the field discontinuity and the average excitation 
via the generalized sheet transition conditions (GSTCs) 
[311351 nans]. 


Zse{y) 


Fsm{y) 


1 Ex {r)\^^Q++ Ej; (r)|^_,,Q- 

l gv(rOUo++g,(r-)|,^o- 

2EAr^\z^o-^-EA^\z^o-' 


( 2 ) 


To promote directive radiation towards 0out "we require 
that the aperture (transmitted) fields approximately fol¬ 
low the suitable plane-wave-like relation 

Ex (fOlz->.o+ ~ ZoutHy 

^ kyloF-^ {^T {kt)} 

(3) 

where Wq {y) is the aperture window (envelope) function 
(yet to be determined) and Zout = 1/^out = ??/ cos^out is 
the TE wave impedance of a plane-wave directed towards 

^out ■ 

In previous work, we have shown that if the wave 
impedance and the real power are continuous across the 
metasurface, these aperture fields can be supported by a 
passive lossless HMS (purely reactive Zge and Ygm) [45] . 
The first condition, local impedance equalization, means 
that the total (incident-l-reflected) fields on the bottom 
facet of the metasurface should exhibit the same wave 
impedance, i.e. Ex{f)\^^Q- = 

achieved by setting the reflection coefficient to a Fresnel- 
like form 


k cos gput - P 
k cos 0out + P ’ 


(4) 
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determining the reflected fields everywhere [Eq. 0 ] , fix¬ 
ing our first degree of freedom. To satisfy the second 
condition, local power conservation, we require that the 
aperture window function follows the magnitude of the 
total (incident-l-reflected) fields at z —>• 0“, namely, 


ITo (z/) = |E, (r-)| 


z—>0“ 






1/2 


(5) 


of oblique-angle radiators is addressed in Appendix [A| . 
In this case, the desirable radiation pattern is symmetric 
with respect to the xz plane, thus it is only natural to use 
a laterally-symmetric excitation; subsequently, we set the 
lateral position of the source to be y' = 0 throughout the 
rest of the paper. With this selection, the even modes 
vanish, and the odd modes follow a cosine profile in the 
lateral dimension. 


B. Optimizing the cavity excitation 


where g * is the autocorrelation of the spectral-domain 
function g (kt) [4^ this determines the transmitted helds 
everywhere [Eq. Q] , fixing our second degree of freedom. 

The absolute value operator in the last equality is of 
utmost significance: it indicates that the transmission 
spectrum of the aperture helds follows, up to a square 
root, the power spectral density (PSD) of {r)\^^Q-, 
and not the spectral content of the incident and rehected 
fields. This is directly related to the balanced (plane- 
wave-like) contribution of the electric and magnetic helds 
to the power how that we stipulated in Eq. (|^, and 
results in a signihcantly favourable plane-wave spectrum, 
as will be discussed in detail in the next subsection. 

Finally, we use these semianalytically predicted helds 
[Eq. Q] and the equivalence principle, manifested by the 
GSTCs [Eq. ([^], to calculate the required HMS surface 
impedance, yielding the desirable purely reactive modu¬ 
lation given by. 


^ se jy) ^ Ysm {y) 
Eout Pout 


-cot 

2 


(t>- {y) - {y) 

2 


( 6 ) 


where cj)± (y) = A {y,z)\^^Q± are the phases of the 
stipulated helds just above and below the metasurface 

Sc¬ 
once the general design procedure is established, ap¬ 
plying it to the conhguration of Fig. which includes an 
electric line source at {y',z') surrounded by PEC walls 
at z = — d, y = EL 12, is straightforward: it reduces to 
hnding the corresponding source spectrum. The latter is 
quantized due to the lateral cavity, and includes multiple 
rehections between the HMS at z = 0 and the PEC at 
z = —d [H]; explicitly, 

^ re-t/3(<^+P) _ ^idid+z') 

n——C5ol Q~3^ty _j_ ^ ^jkty 

(7) 

We refer to the sum of the helds corresponding to the 
n, —n terms in the summation as the held of the nth 
mode of the lateral cavity, where n > 0. 

Although this procedure is applicable for any trans¬ 
mission angle, we restrict ourselves from now on to the 
case of broadside radiation 6*out = 0, where the perfor¬ 
mance of shielded and unshielded FP-LWAs is the most 
problematic due to the splitting condition [2] (design 



One of the key differences between the cavity-excited 
HMS antenna and FP-LWAs is that by harnessing the 
equivalence principle we control the individual contri¬ 
butions of the electric and magnetic helds to the how 
of power, expressed by the lateral distribution of the z- 
component of the Poynting vector on the aperture. More 
specihcally, as mentioned in the previous subsection, the 
resultant (transmitted) aperture helds corresponding to 
Eq. ([^ actually follow the square root of the power pro- 
hle dictated by the cavity mode, and not the prohle of 
the cavity fields. This distinction is very important, as 
the power prohle of a standing wave is always positive, 
whereas the field prohle changes signs along the lateral 
dimensions. Hence, the spectral content of the aperture 
helds, which determines the far-held radiation pattern, is 
fundamentally different. 

To illustrate this point, and to elucidate the guide¬ 
lines for optimizing the cavity excitation, we compare 
the helds formed on the device aperture for a shielded 
FP-LWA, where a standard PRS is used, and for a cavity- 
excited HMS antenna with the same excitation. Figure [^ 
presents the spatial prohle of the tangential electric held, 
its spatial Fourier transform, and corresponding radia¬ 
tion patterns (calculated following [50]), for single mode 
excitation of the n = 1 (solid blue lines), n = 9 (dashed 
red lines), and n = 19 (green solid lines) modes, for an 
aperture length of L = lOA. All plots are normalized to 
their maximum, as the radiation pattern is sensitive to 
the variation of the helds, and not to their magnitude. 

As follows from Eq. 0> the spatial prohle of the 
nth-mode aperture held is proportional to cos(n7ry/L) 
for a standard PRS, but for an HMS designed accord¬ 
ing to Eq. it is proportional to [cos (n7ry/L)| [Fig. 
[^a),(d)]. Except for the lowest-order mode n = 1, 
for which the two functions coincide, the diherence in 
the spatial prohle translates into distinctively diherent 
features in the spectral content [Fig. [^b),(e)]. For 
the nth mode, the transmission spectrum of the HMS 
aperture corresponds to the autocorrelation of the PRS 
aperture spectrum, leading to formation of peaks cen¬ 
tered around the second harmonics {kt = ±2n7r/L) and 
DC {kt =0). As both the right-propagating and left- 
propagating components of the standing wave coherently 
contribute to the DC peak, the latter dominates the 
transmission spectrum, and the radiation patterns cor¬ 
responding to the HMS aperture exhibit highly-directive 
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FIG. 2. Comparison between aperture profiles of cavity-excited PRS, as in shielded FP-LWAs, and cavity-excited HMS, and 
their implication on the radiation patterns. Single mode excitations of the n = 1 (blue solid line), n = 9 (red dashed line), 
and n = 19 (green solid line) modes of an aperture of length L = lOA are compared to the multimode excitation corresponding 
to the HMS antenna presented in Fig. [^with L = lOA, z' = —A, and d = 1.61A (black dash-dotted line). (a),(d) Normalized 
spatial profile of the tangential electric field on the aperture. (b),(e) Normalized spectral content of the aperture field; shaded 
region correspond to the visible part of the spectrum. (c),(f) Normalized radiation patterns. Inset: Close-up of the radiation 
pattern around 6^ = 0. 


radiation towards broadside [Fig. [^f)]. In contrast, the 
PRS aperture devices exhibit symmetric (conical) ra¬ 
diation to angles determined by the dominant spectral 
components of the aperture fields [HIS], i.e. towards 
9 = ± arcsin [n\/ (2L)] [Fig. [^c)]. 

The transverse wavenumber kt = tt/L corresponding 
to the lowest-order mode n = 1 is small enough such 
that the two symmetric beams merge |14j , which enables 
the PRS aperture to radiate a single beam at broad¬ 
side. Indeed, small-aperture shielded FP-LWAs utilize 
this TEio mode to generate broadside radiation. How¬ 
ever, as demonstrated by [50], the aperture efficiency of 
this mode is inherently limited to about 80%, due to the 
non-optimal cosine-shaped aperture illumination of the 
lowest-order mode |55H50] . leading to broadening of the 
main beam [inset of Fig. [^f)]. This highlights a key 
benefit of using an HMS-based antenna, as it is clear 
from Fig. [^f) that we can use high-order mode excita¬ 
tions, which provide a more uniform illumination of the 
aperture, for generating narrow broadside beams with 
enhanced directivities. 

In fact, as the index n of the mode exciting the HMS in¬ 
creases, the autocorrelation of Eq. drives the second 
harmonic peaks outside the visible region of the spec¬ 


trum [shaded region in Fig. ib) ,(e)], funnelling all the 
radiated power to the broadside beam, subsequently in¬ 
creasing the overall directivity. This improvement in ra¬ 
diation properties can be explained using ordinary array 
theory. As seen from Fig. Sd), the peaks of the field 
profile generated by the nth mode on the HMS aperture 
form hot spots of radiating currents separated by a dis¬ 
tance of L/n. The radiation from such an aperture pro¬ 
file would resemble the one of a uniform array with the 
same element separation. As known from established ar¬ 
ray theory, to avoid grating lobes the element separation 
should be smaller than a wavelength HUj. For an aper¬ 
ture length of L = NX, where N is an integer, the hot 
spot separation satisfies this condition for mode indices 
n > N; specifically, for A^ = 10 (Fig. [^, grating lobes 
would not be present in the radiation pattern for mode 
indices n > 10. In agreement with this argument, Fig. 
I^f) shows that for n = 9 grating lobes still exist, while 
for the highest-order mode n = 19 they indeed vanish. 

These observations are summarized in Fig. [^ where 
the radiation characteristics of an HMS aperture of L = 
lOA excited by a single mode are plotted as a function 
of the mode index n (only fast modes kt^n < k are 
considered). For comparison, solid lines denoting the 
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FIG. 3. Radiation characteristics of different lateral cavity 
modes. Half-power beamwidth (blue circles) and 2D directiv¬ 
ity (red asterisks) of an HMS aperture of length L = lOA 
excited by a single mode as a function of the mode index 
n. Solid lines denote the respective radiation characteristics 
of a uniformly excited aperture [30] and dash-dotted lines 
mark the HPBW (blue) and directivity (red) of multimode 
excitation corresponding to the HMS antenna of Fig. [^with 
L — lOA, z' = —A, and d = 1.61A. 


half-power beamwidth (HPBW) and directivity values 
achieved for a uniformly-illuminated aperture of the same 
size are presented as well. Indeed, it can be seen that 
the lowest-order lateral mode exhibits the worst perfor¬ 
mance by far, and the performance improves as the mode 
index increases. While the half-power beamwidth satu¬ 
rates quickly, the directivity values continue to increase 
with n until the point in which grating lobes disappear 
n = N = 10 is crossed; for mode indices n > 10 the radi¬ 
ation characteristics of the HMS aperture are comparable 
with those of the optimal uniformly-excited aperture. 

From an array theory point of view, excitation of the 
highest-order mode is preferable, as the corresponding 
equivalent element separation approaches A/2, implying 
that such aperture profile would be suitable for directing 
the radiation to large oblique angles 0out ^ 0 without 
generating grating lobes [10]. Another reason to pre¬ 
fer excitation of the highest-order mode in the case of 
cavity-excited HMS antennas is that the HMS reflection 
coefficient F (kt) grows larger with kt = wnjL [Eq. @]; 
therefore, the power carried by the highest-order mode 
n = 2N — 1 is best-trapped in the cavity, guaranteeing 
uniform illumination even in the case of very large aper¬ 
tures. 

Nevertheless, generating a single-mode excitation of a 
cavity via a localized source can be very problematic 
[SSI ED. Fortunately, the cavity-excited HMS antenna 
can function very well also with multimode excitation, 
as long as high-order modes dominate the transmission 
spectrum. This is demonstrated by the dot-dashed lines 
in Figs. [^and[^ corresponding to a multimode excita¬ 
tion generated by the configuration depicted in Fig. [T] 
with L = lOA, z' = —A, and d = 1.61A. 

As expected from the expression for the source spec¬ 


trum [Eq. 0]. for a given aperture length L, the field 
just below the aperture due to a line source would be 
a superposition of lateral modes, the weights of which 
are determined by the particular source configuration, 
namely the cavity thickness d and source position z'. The 
multimode transmission spectrum in Fig. 0b) indicates 
that for the chosen parameter values, high-order modes 
{kt —>■ ik) predominantly populate the aperture spec¬ 
trum, however low-order modes {kt —> 0) are present 
as well, to a non-negligible extent (this takes into ac¬ 
count the fact that the transmission coefficient 1 — F (fct) 
is higher for lower-order modes). Considering that the 
far-held angular power distribution S {&) is proportional 
to cos^ 6* |r {kt = A; sin 0)1^, the multimode excitation of 
the PRS aperture results in a radiation pattern resem¬ 
bling the one corresponding to single mode excitation 
of the highest-order mode {n = 19) but with signihcant 
lobes around broadside [Fig. 0c)]; consequently, even if 
a conical beam is desirable, multimode excitation would 
result in signihcant deterioration of the directivity. 

On the other hand, the same multimode excitation 
does not degrade substantially the performance of the 
HMS antenna. The autocorrelated spectrum, relevant to 
the held induced on the HMS aperture, results in merging 
of all spectral components into a sharp DC peak, with the 
high-order grating lobes pushed to the evanescent region 
of the spectrum [Fig. 0e)]. This retains a beamwidth 
comparable with that resulting from a single-mode exci¬ 
tation of the highest-order mode, with only slight deteri¬ 
oration of the directivity due to increased side-lobe level 
[Fig. I^f) and inset]. Continuing the analogy to array 
theory, such multimode excitation introduces slight vari¬ 
ations to the magnitude of the array elements, forming 
an equivalent non-uniform array [10] . The correspond¬ 
ing multimode HPBW and directivity values are marked, 
respectively, by blue and red dash-dotted lines in Fig. 
[^ verifying that indeed, cavity-excited HMS antennas 
achieve near-unity aperture efficiencies with a practical 
multimode excitation; this points out another key advan¬ 
tage of the cavity-excited HMS antenna with respect to 
shielded FP-LWAs. 

With these observations in hand, we are finally ready 
to formulate guidelines for optimizing the cavity excita¬ 
tion for maximal directivity. For a given aperture length 
L = NX, with respect to Eq. Q, we maximize the cou¬ 
pling to the n = 2N — 1 mode (which exhibits the best 
directivity) by tuning the cavity thickness d as to mini¬ 
mize the denominator of the corresponding coupling co¬ 
efficient; equally important, we minimize the coupling to 
the 71 = 1 mode (which exhibits the worst directivity) by 
tuning the source position z' as to minimize the numera¬ 
tor of the corresponding coupling coefficient. To achieve 
these with minimal device thickness we derive the follow¬ 
ing design rules 


A 2N 
2 - 1 






( 8 ) 


Although this is somewhat analogous to the typical de- 
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sign rules for (unshielded) FP-LWAs [13], the key dif¬ 
ference is that for HMS-based antennas we optimize the 
source conhguration regardless of the desirable transmis¬ 
sion angle 0out- This difference is directly related to the 
utilization of the equivalence principle for the design of 
the proposed device, which provides certain decoupling 
between its excitation and radiation spectrum [cf. Fig. 
Sb) ,(e)]. This decoupling becomes very apparent when 
the HMS antenna is designed to radiate towards oblique 
angles 0out 7 ^ 0, in which case the same cavity excitation 
yields optimal directivity as well (See Appendix [A|) . 

Two important comments are relevant when consid¬ 
ering these design rules. First, even though following 
Eq. ([^ maximizes the coupling coefficient of the highest- 
order mode and minimizes the coupling coefficient of the 
lowest-order mode, it does not prohibit coupling to other 
modes. The particular superposition of lateral modes ex¬ 
hibits a tradeoff between bemawidth and side-lobe level 
(as for non-uniform arrays [lOjl. Thus, final optimiza¬ 
tion of the cavity illumination profile is achieved by fine- 
tuning the source position z' for the cavity thickness d de¬ 
rived in Eq. Q , with the aid of the efficient semianalyti- 
cal formulas. In fact, the source position z' is another de¬ 
gree of freedom that can be used to optimize the radiation 
pattern for other desirable performance features, such as 
minimal side-lobe level; this feature is further discussed 
and demonstrated in Appendix]^ Second, although the 
optimal device thickness increases with increasing aper¬ 
ture length, the increase is sublinear, following an asymp¬ 
totic square-root proportion factor. Therefore, applying 
the proposed concept to very large apertures would still 
result in a relatively compact device, while efficiently uti¬ 
lizing the aperture for producing highly-directive pencil 
beams. 


III. RESULTS AND DISCUSSION 

We follow the design procedure and the considerations 
discussed in Section |n| to design cavity-excited HMS an¬ 
tennas for broadside radiation with different aperture 
lengths: L = lOA, L = 14A, and L = 25A. The cavity 
thickness was determined via Eq. Q to be d = 1.61 A, 
d = 1.89A, and d = 2.50A, respectively; the source 
position was set to z' = —l.OOA, z' = —1.33A, and 
z' = — 1.94A, respectively, exhibiting maximal directiv¬ 
ity. 

The required electric surface impedance and magnetic 
surface admittance modulations [Eq. (|^] are imple¬ 
mented using the ’’spider” unit cells depicted in Fig. 

At the design frequency / = 20GHz (A « 15mm), the 
unit cell transverse dimensions are A/10 x A/10 and the 
longitudinal thickness is 52mil « A/12. Each unit cell 
consists of 3 layers of metal traces defined on two bonded 
laminates of high-dielectric-constant substrate (See Ap¬ 
pendix]^. The two (identical) external layers provide 
the magnetic response of the unit cell, corresponding to 
the magnetic surface susceptance B^rn = which 



FIG. 4. Physical configuration of the ’’spider” unit cells used 
for implementing the HMS at a frequency of / = 20GHz (A « 
15mm). The electric response is controlled by the capacitor 
width H4 of the electric dipole, while the magnetic response 
is determined by the magnetic dipole arm length (See 
Appendix [C] for detailed description). 

is tuned by modifying the arm length Lm (affects mag¬ 
netic currents induced by tangential magnetic fields Hy). 
Analogously, the middle layer is responsible to the elec¬ 
tric response of the meta-atom, corresponding to the elec¬ 
tric surface reactance Xge = 9{-^se}, which is tuned by 
modifying the capacitor width We (affects electric cur¬ 
rents induced by tangential electric fields Ex). By con¬ 
trolling the lengths of and Wg, the spider unit cells 
can be designed to exhibit Huygens source behaviour, 
with balanced electric and magnetic responses ranging 
from BsmT] = Xseh = -3.1 to Bsmll = Xseh = 0-9 
(See Appendix [C| . 

Figure presents the design specifications, field dis¬ 
tributions, and radiation patterns for the three cavity- 
excited HMS antennas; Table [l| summarizes the antenna 
performance parameters (for reference, parameters for 
uniformly-excited apertures |3D| are also included). The 
semianalytical predictions [Eq. 01 are compared to full- 
wave simulations conducted with commercially-available 
finite-element solver (ANSYS HFSS), where the HMS 
was implemented using the aforementioned spider cells 
(See Appendix 0. As demonstrated by Fig. [ll;a)-(c), 
the realized unit cells are capable of reproducing the 
required surface impedance modulation, except maybe 
around large values of Bsmil = A^e/?!; however, it has 
been shown that such discrepancies usually have little 
effect on the performance of Huygens’ metasurfaces [52j . 










TABLE 1. Radiation characteristics of cavity-excited HMS antennas (corresponding to Fig. [^. 



L = lOA (d = 

1 . 6 IA, 12 ' 

= l.OOA) 

L = 14A (d = 

1.89A, |z' 

= 1.33A) 

L = 25A (d = 

2.50A, |z' 

= 1.94A) 


Full-wave Semianlytical Uniform 

Full-wave Semianlytical Uniform 

Full-wave Semianlytical Uniform 

HPBW 

5.11° 

5.13° 

5.08° 

3.83° 

3.64° 

3.63° 

2.09° 

2.13° 

2.03° 

Directivity (2D) [dBi] 

17.42 

17.84 

17.98 

18.79 

19.15 

19.44 

21.33 

21.75 

21.96 

First Side-Lobe 

00 

0 

0 

cx 5 

8.2° 

6.4° 

6.3° 

5.9° 

0 

CO 

0 

CO 

3.3° 

Side-Lobe Level [dB] 

- 10.4 

-12.0 

-13.5 

-11.1 

-10.4 

-13.5 

-14.6 

-14.0 

-13.5 
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FIG. 5. Performance of cavity-excited HMS antennas with aperture lengths of L = lOA, L = 14A, and L = 25A. (a)-(c) HMS 
design specifications Xse (y) /7 = Bam iy)v (black solid line) derived from Eq. Q, and the realized electric surface reactance 
(blue circles) and magnetic surface susceptance (red circles) using the spider unit cells, (d)-(f) Radiation patterns produced by 
semianalytical formalism (blue dashed line) and full-wave simulations (red solid line), (g)-(i) Field distribution {y,z)}\ 

produced by full-wave simulations, (j)-(l) Semianalytical prediction of \^{Ea: {y,z)}\ |45| . 
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The results in Fig. [^and Table[T]indicate that the fields 
and radiation properties predicted by the semianalyical 
formalism are in excellent agreement with the full-wave 
simulations for a wide range of aperture lengths. The 
utilization of realistic (lossy) models for the conductors 
and dielectrics in the simulated device, as well as other 
deviations from the assumptions of the design procedure 
(Appendix [^, result in some discrepancies between the 
full-wave simulations and predicted performance; how¬ 
ever, these mostly affect radiation to large angles [Fig. 
[^d)-(f)]. While this contributes to a minor quantitative 
difference in the directivity, the properties of the main 
beam and the side lobes follow accurately the semiana- 
lytical results (Table |T]), indicating that the theory can 
reliably predict the dominant contributions to the radia¬ 
tion pattern, as discussed in reference to Fig. 

For all cases considered, the excitation of the highest- 
order lateral mode is clearly visible [Fig. [^g)-(l)], leading 
to beamwidth and (2D) directivity values D compara¬ 
ble with the ones achieved by uniform excitation of the 
aperture (Table |T|. In particular, the simulated radiation 
patterns yield aperture efficiencies ?7apr = D/ (27rL/A) of 
88%, 86%, and 87% for the L = lOA, L — 14A, and 
L = 25A devices, respectively, retained even when the 
aperture length is very large. In terms of half-power 
beamwidth, often taken as a measure for effective ex¬ 
ploitation of the aperture m, the device performance is 
even closer to that of a unifromly-excited aperture, with 
pencil beam HPBWs reaching 99%, 95%, and 97% of the 
optimal beamwidth, for L = lOA, L — I4A, and L = 25A 
devices, respectively. It should be stressed that even 
though optimized TEio shielded FP-LWAs can reach 
aperture efficiencies of 80%, their HPBWs are limited 
to about 75% of the optimal beamwidth [3D]; more im¬ 
portantly, their PRS-based design requires single-mode 
excitation to achieve this performance, thus preventing 
practical realization of large-aperture devices. 

We would like to emphasize that this near-optimal 
aperture utilization is achieved while using realistic mod¬ 
els for the substances comprising the metasurface, as well 
as geometrical dimensions compatible with practical fab¬ 
rication techniques. It is also very clear from the results 
that there is no apparent degradation of the structure 
performance for increasing aperture length, even for the 
very large values considered. Therefore, it would be rea¬ 
sonable to conclude that the concept introduced in this 
paper can be applied to design cavity-excited HMS an¬ 
tennas with arbitrarily-large apertures which will exhibit 
near-unity aperture efficiency. 


IV. CONCLUSION 

We have introduced a novel design for low-profile 
single-fed antennas exhibiting beamwidth and directivity 
values comparable with uniformly-excited apertures. To 
that end, we harness the equivalence principle to devise a 
cavity-excited Huygens’ metasurface, setting the source 


configuration, HMS reflection coefficient, and aperture 
fields such that (1) the highest-order lateral mode is pre¬ 
dominantly excited, which guarantees that the aperture 
is well-illuminated; (2) the aperture fields follow the in¬ 
cident power profile and not the incident field profile, 
which forms an array-like aperture profile with favourable 
transmission spectrum; and (3) the power flow and wave 
impedance are continuous across the metasurface, which 
ensures the design can be implemented by a passive and 
lossless HMS. The possibility to control the field disconti¬ 
nuities using the electric and magnetic currents induced 
on the HMS allows us to optimize separately the cav¬ 
ity excitation and the radiated fields, thus overcoming 
the fundamental tradeoff existing in FP-LWAs between 
aperture efficiency and edge-taper losses. 

It should be emphasized that the general design pro¬ 
cedure formulated and demonstrated herein facilitates 
further optimization of such devices for various appli¬ 
cations. The extensive freedom one has in choosing the 
source configuration, combined with the efficient semi- 
analytical approach, allows explorations of other excita¬ 
tion sources, e.g. with different orientations and current 
distributions, to tailor the aperture fields for one’s re¬ 
quirements, e.g. minimal side-lobe level or radiation to¬ 
wards an oblique angle (See Appendices [A| and [b|) ; once 
the source spectrum is characterized, the rest of the pro¬ 
cedure is straightforward, and the fields and radiation 
patterns are readily predicted. In particular, the for¬ 
malism can be readily applied to design cavities based 
on dielectric Bragg reflectors instead of metallic mirrors, 
more suitable for terahertz and optical devices. In view 
of the recent demonstrations of terahertz and optical 
(Huygens’) metasurfaces |5D11331 [53H55] . this would al¬ 
low realization of compact and efficient pencil beam ra¬ 
diators across the electromagnetic spectrum, extending 
the range of applications even further, e.g. for efficient 
outcoupling of single-photon emission |59j or terahertz 
elementary sources [60] . 


Appendix A: Cavity-excited HMS antennas 
radiating at oblique angle 


To design cavity-excited HMS antennas radiating at 
oblique angle we follow the same procedure outlined in 
Section |IT| for broadside radiators, with the desirable 
^out 0 inserted into Eqs. ^ and Q; note that now 
Zout = l/Lbut = Jy/cos^out ^ V- The same argumen¬ 
tations for optimizing the cavity excitation holds (Sub¬ 
section IIBj ), where the HMS transmission spectra plot¬ 
ted in Fig. [^e) merely shifted by kt = fcsin^out- This 
is essentially identical to the design of antenna arrays 
radiating at an angle 0outj where the magnitude of the 
array element currents are usually determined indepen¬ 
dently of 0out, and the direction of the main beam is 
manifested via phase-shifts imposed between adjacent el¬ 
ements HDj. Therefore, the optimal configuration still 
follows Eq. (|^, i.e. we use the same cavity excita- 
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tion regardless of 0out ■ Once the required electric surface 
impedance and magnetic surface susceptance are spec¬ 
ified, the corresponding spider unit-cell dimensions are 
retrieved from a lookup table constructed as described in 
Appendix with Zout = 1/i^out = ?7/cos6>out- 

We use this outlined procedure to design cavity-excited 
HMS antennas radiating towards 0out = 30°. The de¬ 
sign specifications and radiation characteristics are pre¬ 
sented in Fig. and Table [IT] for devices with aperture 
lengths of L = lOA and L = 14A, comparing the re¬ 
sults of full-wave simulations, semianalytical predictions, 
and uniformly-excited apertures (note that for the lat¬ 
ter, HPBW and directivity are factored by cos 0out with 
respect to their values for broadside radiation [ID]). In 
consistency with the aforementioned decoupling between 
the excitation and main-beam angle, the optimal cavity 
configuration for maximal directivity is identical to the 
one derived for broadside radiation (compare Tables □ 
and 0- Indeed, comparison of Fig. |^e)-(h) with Fig. 
indicates that the field profile inside the cavity is practi¬ 
cally independent of 0out- 

To match better the range of required surface 
impedance values [Fig. [^a)-(b)] to the one achievable by 
our spider unit-cells (Appendix [Cj) , we introduce a con¬ 
stant phase-shift ^out to the aperture fields [DD] HSj . This 
adds a constant ■fout/2 to the argument of the cotangent 
surface impedance modulation [Eq. (|^], which varies 
the distribution of required surface impedance values, 
but does not affect the radiation pattern. Simulation 
methodology follows Appendix [^ with the difference 
that the PMC in the xz plane cannot be used, as the 
HMS is not symmetric in this case [Fig. [^a)-(b)]. This 
doubles the volume of the simulation domain, resulting in 
insufficient convergence of the computations for L = 25A 
devices. Hence, only results for L = lOA and L = 14A 
are presented. 

Table [H| and Fig. [^indicate that the excellent agree¬ 
ment between the semianalytical predictions and full- 
wave simulations is retained for devices radiating at 
oblique angles. For the L = lOA and L = 14A antennas, 
respectively, very high aperture efficiencies of Tyapr = 83% 
and r?apr = 85% are recorded. The HPBWs are within 
94% and 101% of the optimum corresponding to uniform 
illumination, comparable to the characteristics achieved 
for broadside-radiating cavity-excited HMS antennas. 



FIG. 6. Performance of cavity-excited HMS antennas radiat¬ 
ing towards 6^out = 30°, with aperture lengths of L = lOA and 
L = 14A. (a)-(b) Specified (blue dashed line) and realized 
(blue open circles) electric surface reactance, and specified 
(red solid line) and realized (red open circles) magnetic sur¬ 
face susceptance using the spider unit-cells, (c)-(d) Radiation 
patterns produced by semianalytical formalism (blue dashed 
line) and full-wave simulations (red solid line), (e)-(f) Field 
distribution \^{Ex {y, z)}\ produced by full-wave simulations. 
(g)-(h) Semianalytical prediction of |K {Aa, (y,z)}\ |45| . 


Appendix B: Cavity-excited HMS antennas with 
reduced side-lobe level 


As denoted in Subsection HB[ for an optimized cav¬ 
ity thickness d, i.e. one that maximizes coupling to the 
highest-order mode, one may utilize the source position 
z' as an additional degree of freedom to optimize the an¬ 
tenna radiation characteristics to achieve desirable per¬ 
formance. In Sections |H| and [Hi] we utilized this degree 
of freedom to suppress the coupling to the lowest-order 
mode, in order to facilitate the highest possible directiv¬ 


ity. This, in consistency with array theory m , comes 
at the expense of side-lobe level. To demonstrate the 
possibility to utilize the efficient semianalytical formu¬ 
lation to devise a cavity excitation which reduces the 
side-lobe level, we sweep the source position z' for the 
optimal d = 1.61 A corresponding to an aperture length 
of L = lOA, and evaluate the directivity and side-lobe 
level associated with each z'. As we only require exami¬ 
nation of the radiation pattern properties, it is sufficient 
to compute the radiated fields, which can be achieved 
analytically by asymptotic evaluation of Eq. 0 in con- 
junction with Eq. Q (See, e.g. Appendix C of |i5]l. 
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TABLE II. Radiation characteristics of cavity-excited HMS antennas radiating towards 0out = 30° 
(corresponding to Fig. [^. 



{d = 1.61A 

L = lOA 

\z'\ = 1.00A,^out = ^tt) 

(d= 1.89A 

L = 14A 

\z'\ = 1.33A,^out = ^rr) 


Full-wave Semianlytica 

Uniform 

Full-wave Semianlytica 

Uniform 

Main beam 

30.2° 

30.0° 

30.0° 

30.2° 

30.0° 

30.0° 

HPBW 

6 .20° 

5.90° 

5.86° 

H-' 

to 

0 

4.20° 

4.20° 

Directivity (2D) [dBi] 

16.54 

17.23 

17.36 

18.12 

18.49 

18.82 

Side lobe 

20.55° 

20.49° 

20.49° 

23.32° 

22.96° 

23.13° 

Side-Lobe #1 Level [dB]^ 

-13.7 

-11.4 

-13.5 

-9.0 

-9.7 

-13.5 

Side lobe #2^ 

41.04° 

40.54° 

40.54° 

37.42° 

37.59° 

37.38° 

Side-Lobe #2 Level [dB]^ 

-12.0 

-12.3 

-13.5 

-10.2 

-10.4 

-13.5 


^ Side lobes #1 and #1 refer, respectively, to the first side lobes at angles lower and higher than 
the main beam angle. 


The results of the parametric sweep are presented in 
Fig. 0 with indication of the maximal directivity con¬ 
figuration {\z'\ = lA) and a conhguration with reduced 
side-lobe level [\z'\ = 1.25A). The design specifications 
and radiation characteristics of the cavity-excited HMS 
antenna corresponding to the latter are presented in Fig. 
0 and Table |III| As expected from array theory, the re¬ 
duction of side-lobe level (by about lOdB according to 
full-wave simulations) results in a broadening of the main 
beam and a reduction of the overall aperture efficiency 
(although still quite high value of rjapr = 81% is recorded 
by full-wave simulations). This demonstrates the versa¬ 
tility of our approach, with which a range of radiation 
pattern properties can be achieved. 



FIG. 7. Utilizing the source position as an additional degree of 
freedom for achieving desirable radiation characteristics. Di¬ 
rectivity (red solid line) and maximal side-lobe level (dashed 
blue line) of a cavity-excited HMS antenna with L — lOA and 
d — 1.61A are presented as a function of source position \z'\. 
The source position maximizing directivity {\z'\ = lA, corre¬ 
sponding to Fig. 0 is denoted by a red dotted line, and the 
source position for reduced side-lobe level is denoted W a blue 
dash-dotted line {\z'\ = 1.25A, corresponding to Fig. 0. 


TABLE III. Radiation characteristics of a cavity-excited 
HMS antenna with reduced side-lobe level (corresponding 
to Fig. [^. 



L = lOA (d = 

1.61A, ja' 

= 1.25A) 


Full-wave Semianlytical Uniform 

HPBW 

6.13° 

5.89° 

5.08° 

Directivity (2D) [dBi] 

17.04 

17.62 

17.98 

First Side lobe 

9.62° 

9.21° 

8 .2° 

Side-Lobe Level [dB] 

-20.7 

-20.5 

-13.5 


Appendix C: Spider unit-cell modelling 

The spider unit cells depicted in Fig. 0were defined 
in ANSYS Electromagnetic Suite 15.0 (HFSS 2014) with 
two 25mil-thick (« 0.64mm) Rogers RT/duroid 6010LM 
laminates (green boxes in Fig. 0 bonded by 2mil- 
thick (« 51/im) Rogers 2929 bondply (white box in Fig. 
0. The electromagnetic properties of these products at 
20GHz, e.g. permittivity tensor and dielectric loss tan¬ 
gent, as were provided to us by Rogers Corporation, have 
been inserted to the model. Specifically, a uniaxial per¬ 
mittivity tensor with e^x = ^yy = ^3.3eo,ezz = lO.Sleo 
and loss tangent of tan 5 = 0.0023 were considered for 
Rogers RT/duroid 6010LM laminates, while an isotropic 
permittivity of e = 2.94eo and loss tangent tani5 = 0.003 
were considered for Rogers 2929 bondply. The copper 
traces corresponded to 1/2 oz. cladding, featuring a thick¬ 
ness of 18/rm; the standard value of cr = 58 x 10®S/m 
bulk conductivity was used in the model. To comply 
with standard PCB manufacturing processes, all copper 
traces were 3mil (« 76/im) wide, and a minimal distance 
of 3mil was kept between adjacent traces (within the cell 
or between adjacent cells). This implies that the hxed 
gaps between the capacitor traces (along the x axis) of 
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FIG. 8. Performance of a cavity-excited HMS antenna with 
reduced side-lobe level {L = 10A), (a) HMS design specifi¬ 
cations Xse (y) jf] = Bsra {v )(black solid line) derived from 
Eq. and the realized electric surface reactance (blue cir¬ 
cles) and magnetic surface susceptance (red circles) using the 
spider unit cells, (b) Radiation patterns produced by semi- 
analytical formalism (blue dashed line) and full-wave simu¬ 
lations (red solid line), (c) Field distribution |5R (?/, 2)}| 

produced by full-wave simulations, (d) Semianalytical predic¬ 
tion ot\'i!i{E^{y,z)}\ [4^. 


the electric dipole in the middle layer, as well as between 
the two arms (along the y axis) of the magnetic dipole in 
the top and bottom layer (Fig. |^, were fixed to a value 
of Dg = 3mil (« 76/im); the distance from the arm edge 
to the edge of the unit cell was fixed to Dg/2 = 1.5mil 
(« 38/rm). 

Unit cells with different values of magnetic dipole arm 
length Lm and electric dipole capacitor width We were 
simulated using periodic boundary conditions; HFSS Flo- 
quet ports were placed at z = ±A and used to character¬ 
ize the scattering of a normally-incident plane wave off 
the periodic structure (the interface between the bond- 


ply and the bottom laminate was defined as the z = 0 
plane). For each combination of and We, the corre¬ 
sponding magnetic surface susceptance and electric 
surface reactance Xgg were extracted from the simulated 
impedance matrix of this two-port configuration, follow¬ 
ing the derivation in [61j . 

The magnetic response Bgrn was found to be propor¬ 
tional to the magnetic dipole arm length Lm, with al¬ 
most no dependency in Wg [10]. Thus, to create an 
adequate lookup table for implementing HMSs radiat¬ 
ing towards 6>out, we varied Lm by constant increments, 
and for a given Lm, plotted Bgrn/Yont and Xge/Zout 
as a function of Wg- The value of Wg for which 
the two curves intersected corresponded to a balanced- 
impedance point (ZgelZoe^t = Ysm/Yout), where the unit 
cell acts as a Huygens source, and thus suitable for im¬ 
plementing our metasurface. A lookup table composed 
of {Bsrn, Xgg) pairs and the corresponding unit cell ge¬ 
ometries {Lm,Wg) was constructed, and refined through 
interpolation. The interpolated unit cell geometries were 
eventually simulated again, to verify the interpolation ac¬ 
curacy and finalize the lookup table entries, as presented 
in Fig. [^ Finally, for a given HMS with prescribed 
surface impedance modulation {Bg^ {y) ,Xgg (y)), a cor¬ 
responding structure could be defined in HFSS using the 
unit cells {Lm (y), Wg {y)) found via the lookup table in 
terms of least-squares-error. 



FIG. 9. Graphic representation of spider unit-cell lookup ta¬ 
bles corresponding to HMSs radiating towards Oout = 0° and 
^out = 30°. Gapacitor width values Wg required for achiev¬ 
ing balanced electric and magnetic responses Bsm/Tout = 
Xse/Zout are presented as a function of the magnetic dipole 
arm length Lm for 0out = 0° (blue open squares) and 
Oout = 30° (red open triangles) radiators, as obtained by 
finite-element simulations. The corresponding Bsm/Yout = 
XsgjZogA values are denoted using a blue dashed line for 
fi'out = 0° {Zout = l/Uout = y) and using a red solid line 
for 0out = 30° {Zout = 1/Uout = v/ cos 30°). 
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Appendix D: Antenna full-wave simulations 


To verify our semianalytical design via full-wave sim¬ 
ulations, each of the cavity-excited HMS antennas de¬ 
signed in this paper was defined in HFSS using a sin¬ 
gle strip of unit cells implementing the metasurface, oc¬ 
cupying the region |a:| < A/20, \y\ < Lj^ {L being 
the aperture length of the antenna), and —0.64mm < 
z < 0.69mm (in correspondence to the laminate and 
bondply thicknesses). The simulation domain included 
kl < A/20, \y\ < L/2 + 2.5A, and -d < z < lOA (d 
being the cavity thickness), where PEC boundary con¬ 
ditions were applied to the x = ±A/20 planes to form 
the equivalence of a 2D scenario. PEC boundary condi¬ 
tions were also applied to the z = —d plane, and to two 
18^m-thick side-walls at y = forming the cavity. 

The line-source excitation was modelled by a A/20-wide 
lA current sheet at z = z/ with the current aligned with 
the X axis. Radiation boundary conditions were applied 
to the rest of the simulation space boundaries, namely 
z = lOA, and y = ± (L/2-|-2.5A), allowing proper nu¬ 
merical evaluation of the fields surrounding the antenna. 

To reduce the computational effort required to solve 
this configuration, we utilized the symmetries of our 
TE scenario. Specifically, we placed a perfect-magnetic- 
conductor (PMC) symmetry boundary conditions at the 
xz plane, and a PEC symmetry boundary conditions at 
the yz plane (the PMC symmetry boundary conditions is 
only applicable for broadside radiators, cf. Appendix [A|. 
We also noticed that adding a thin layer (Imil « 25/rm) 
of copper between the electric dipole edges and the PEC 
parallel-plates at a; = ±A/20 enhanced the convergence of 
the simulation results. With that minor modification, all 
of the simulated antennas converged within less than 40 
iterations (maximum refinement of 10% per pass), where 
the stop conditions was 3 consecutive iterations in which 
AEnergy < 0.03. 


Appendix E: Design procedure assumptions 


Several assumptions made during the derivation of the 
HMS design formulas contribute to discrepancies between 
predicted and actual performance of the presented an¬ 
tennas. First, the predicted fields are derived assuming 
the HMS is capable of implementing continuous surface 
impedance boundary conditions, with unbound surface 
impedance values; nevertheless, the physical implemen¬ 
tation requires discretization of the continuous modula¬ 
tion into unit-cells, and the range of achievable surface 
impedance values is limited (See Fig. [^. Second, the 
HMS is assumed to be passive and lossless, however real¬ 
istic conductors and dielectrics, used for the implementa¬ 
tion of the devices in ANSYS HFSS, include unavoidable 
losses. Third, to facilitate the plane-wave-like relation 
between the transmitted fields on the aperture [Eq. ([^], 
while still guaranteeing they obey Maxwell’s equations, 


we have used the approximation 




' {^0Tikt)}\ « |j-i {i-^Tih) [1 ±r(fc*)]} 


which is satisfied when 


Hy) = 






< 1 , 


(El) 

(E2) 


which is a refinement of the slowly-varying envelope 
(SVE) approximation utilized in [45) . This approxima¬ 
tion is self-consistent with our design scheme, as when 
the transmitted fields are directive towards 0out as desir¬ 
able, the dominant components of the transmission spec¬ 
trum T (kt) are in the vicinity of kt = fcsin^out, where 
the reflection coefficient F (kt) completely vanish (the nu¬ 
merator of Eq. (E2) vanishes). 

Interestingly, the impacts of these three assumptions 
can be assessed by reviewing the predicted and simulated 
power flow across the metasurface. The z-directed power 
profiles below (blue) and above (red) the HMS as pre¬ 
dicted by the semaianalytical formalism (open circles and 
squares, respectively) and as extracted from full-wave 
simulations (dashed and solid lines, respectively) are pre¬ 
sented in Fig. for the three antennas reported in Sec¬ 
tion [l^ The fact that the general trend and quantitative 
data of the semianalytical and simulated results compare 
well (note that the profiles are plotted using a common 


77 (/o/A)^ unit scale), indicates that the first assumption 
is valid (this is also supported by [S2]). The semianalyt¬ 
ical predictions made based on a homogenized continu¬ 
ous surface impedance boundary conditions mostly agree 
with the simulation data recorded below and above the 
metasurface, where effective medium theory predicts dis¬ 
cretization effects to be negligible [13117]. 

Violations of the second assumption, regarding the 
lossless nature of the HMS, would manifest themselves 
as differences between the simulated power profile below 
and above the metasurface, which must originate in dissi¬ 
pation in the unit-cell conductors and dielectrics. On the 
other hand, violations of the third assumption, related to 
the SVE approximation, would manifest themselves as 
differences in the semianalytically predicted power pro¬ 
file below and above the metasurface, as they correspond 
to violations of local power conservation jUj . 

Although local deviations from these two assumptions 
are found to be rather small, they contribute to a non- 
negligible reduction of the total power flow across the 
metasurface (integrated over the aperture length). The 
values denoted in the legends of Fig. [^indicate that ac¬ 
cording to full-wave simulations about 15% of the power 
available below the HMS is dissipated in the lossy con¬ 
ductors and dielectrics, while the semianalytically predic¬ 
tions reveal about 10% discrepancy between the power 
below and above the metasurface. 

While these relative deviations can be considered small 
albeit non-negligible, it seems that they actually bal¬ 
ance each other. The theoretical derivation assumes 
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Power profile (L=14(X]) 



Power profile (L=25t>.]) 


-P*t„,(Full-wave) = 0.596769 [t|*I^/ 1 ] 

□ P* “'(Semlanalyllcal|=0.7464S2 [ii‘lj/il 

- P; (Full-wave) = 0.72292 [i,*)'/!] 

o P’“'(Semianalyllcal)=0.8212S7 [n'ljft.] 



FIG. 10. Power profiles on the cavity-excited HMS antenna apertures. The real part of the 2 -component of the Poynting vector 
= 5R{5'z [y,z —>■ 0^)} above and below and the metasurface is presented for the antennas reported in Table I and Fig. 
5, with aperture lengths (a) L — lOA (b) L — 14A, and (c) L = 25A. Full-wave simulation results of (red solid line) and 
Pp (blue dashed line) were evaluated at 2 = A/10 and 2 = —A/10, respectively. Semianalytical predictions of P/" (red open 
squares) and Pp (blue open circles) are presented as well. The total power P/*^tot 5 calculated by integrating the power profile 
along the aperture is indicated in the legend. 


and prescribes a lossless HMS, but the minor viola¬ 
tions of the SVE approximation contribute to predicted 
(Maxwellian) fields which must be supported by small 
losses. On the other hand, the implemented HMS does 
include realistically unavoidable losses, which turn out to 


dissipate a comparable amount of power. We hypothe¬ 
size that this balance allows overcoming the minor devi¬ 
ations from the theoretical assumptions, facilitating the 
very good agreement between predicted and simulated 
results reported herein. 
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